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We demonstrate that the transport characteristics of deep optical lattices with one or multiple 
off-resonant external energy offsets can be greatly-enhanced by modulating the lattice depth in 
an exotic way. We derive effective stationary models for our proposed modulation schemes in the 
strongly interacting limit, where only one particle can occupy any given site. Afterwards we discuss 
the modifications necessary to recover transport when more than one particle may occupy the lattice 
sites. For the specific five-site lattices discussed, we numerically predict transport gains for ranging 
from 4.7 x 10 6 to 9.8 x 10 8 . 
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I. INTRODUCTION 

Ultracold atoms in optical lattice potentials are very rich and diverse physical systems. Their coherent dynamics 
make them candidates for quantum logic gates mm- It has also been proposed that atomtronics components [3H5] 
can be realized in Bose-Hubbard lattices. The typical scheme that might generate such components would involve 
customized lattices that exhibit certain transport characteristics. Holographic mask techniques [Bj make it possible to 
realize such systems in optical lattices, while the current ability to image single atoms makes signal detection feasible. 
Optical lattices have also proven to be excellent quantum simulators: not only has the mott-insulator to superfluid 
phase transition been realized [7] , but lattice shaking ISIS], a modulation technique that involves periodically displacing 
the location of the wells, has recently been used to realize the Haldane model m as well as the Meissner effect El 
in these systems. This lattice shaking technique has also been used to realize photon-assisted tunneling [12HT4] , in 
which the modulation is used to facilitate transport in otherwise frustrated optical lattice systems. 

An alternative modulation technique known as amplitude modulation involves varying the depth of the optical lattice 
V ( t ) sinusoidally. This method, referred to as AC-shaking, is also capable of creating synthetic gauge fields in optical 
lattice systems m- This technique was originally developed as a spectroscopic tool to study optical lattices em¬ 
it was then proposed as a means to populate higher Bloch bands of the optical lattice fT8h20] . More recently, multiple 
frequency or polychromatic modulations were used to stimulate two-photon transitions that coherently promoted a 
significant number of atoms from the ground to the fourth excited Bloch band m- Photon-assisted tunneling has 
also been achieved in systems using amplitude modulation [22l 123] , 

This article is concerned with optimizing the transport characteristics of irregular optical lattices with external 
energy offsets that would otherwise greatly reduce the transport within the lattice. The work presented can be viewed 
as a generalization of photon assisted tunneling with amplitude modulation. Here we study optical lattices where 
V(t) is not necessarily varied perturbatively or sinusoidally, and there are potentially several unique energy offsets 
within the lattice structure. The immediate systems that we have in mind are customized optical lattices, produced 

(a) 



FIG. 1. An illustration of a one-dimensional optical lattice whose end sites are driven by ultracold atomic reservoirs. In the 
strongly-interacting case where the reservoirs are conditioned to put only one atom into the system, the energetically flat lattice 
depicted in (a) yields an optimal transport response, whereas the lattice depicted in (b) yields virtually zero transport. 

possibly by holographic masks and/or other techniques that can generate such lattices [24], that are designed to 
have exotic transport characteristics. For instance, the development of synchronous atomtronics logic elements, may 
require transport across coupled lattice sites with energy mismatches that are large compared to the tunneling rates. 
However, the results we present are very general, and they can be applied to a variety of systems: although we analyze 
the one-dimensional open bosonic case in this work, a variation of the presented schemes are valid for either bosons 
or fermions, trapped in open or closed quantum systems comprised in one-, two-, or three-dimensional optical lattices 
or possibly atom-chip experiments. They might find applications in an experiment where the control of a coherent 
atomic signal is desired. Finally, the modulation schemes might be useful as a spectroscopic tool: scanning through 
modulation frequencies of an arbitrary lattice may lead to spikes in the transport response, which could be used 
to infer properties of the structure of the lattice itself. Such a tool may be useful for studies of disordered optical 
lattices [25 27 . 

We demonstrate that the transport response of irregular, deep, lowest band lattices can be significantly enhanced by 
a polychromatic modulation of the tunneling rates of the system. As discussed below, this requires a unique periodic 
modulation of the lattice depth itself. In a regime where the external energy differences are large compared to the 
intrasystem tunneling, we demonstrate that this modulation enhances transport by mapping the formerly stationary 
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irregular system Hamiltonian that prohibits transport onto an effective stationary Hamiltonian that accurately predicts 
the optimized transport response. 

We examine the transport properties of these systems by driving their end sites with reservoirs of neutral, ultracold 
atoms. When a chemical potential difference exists between the reservoirs, atomic transport (or current) may be 
induced across the system. For example, when the reservoirs and the lattice are conditioned so that each system site 
can support at most one atom, the flat energy landscape depicted in Fig. HJa) allows the atoms to freely explore the 
entire lattice. This configuration yields the optimal transport response. If the landscape has one or more external 
energy differences, as in Fig. [ljb), the current can be greatly suppressed. The schemes presented in this work map a 
system such as the one illustrated in Fig. [ljb) onto an effective lattice of the type in Fig. [lja). 

The conditions for transport optimization turn out to depend on the maximum site occupancy, which is determined 
by the value of the driving chemical potential. After examining the case where a maximum of one atom may exist on 
any given site, we follow with a discussion of current recovery for situations involving higher site occupancies. 

The structure of this article goes as follows: in Sec. [IT] we introduce the time-dependent system Hamiltonian, discuss 
our open quantum system theory and define the observable that we refer to as ‘transport current’. In Sec. m we 
analyze the single site occupancy case. A unitary transformation is used to demonstrate the formal equivalence between 
stationary flat lattices and arbitrary lattices whose tunneling rates undergo a specific complex phase modulation. We 
then adapt this reasoning to propose an amplitude modulation scheme that leads to large gains in an irregular 
lattice’s transport response. Employment of a secular approximation allows us to derive effective stationary models 
that accurately reproduce the current dynamics of these complicated offset modulated systems. Afterwards, numerical 
simulations are provided to confirm and characterize the recovered transport response. In Sec.[lV]we discuss transport 
recovery for the case where up to two atoms can occupy any given lattice site. Extensions to higher site occupancies 
is straight-forward. 


II. THEORY 
A. System Hamiltonian 

The systems we consider in this article are ultracold bosonic atoms trapped in an irregular one-dimensional optical 
lattices whose lattice depth V(t) can be modulated in time. The lattice potential is assumed to be Vi at (x,t) = 
V(t) sin 2 ( 7 rx/ai) + V ext (x ), where ai is the lattice spacing and V ex t, such that \V ext (x)\ <C V(t) for all x and £, is the 
external potential that generates the energy shifts in the lattice structure. 

Working in the lowest Bloch band and in the tight-binding regime, the stationary system dynamics are well-described 
by the Bose-Hubbard model [28| . If the modulation timescale is slow compared to the system dynamics, then the 
time evolution of the quantum states can be treated adiabatically EDI- Under these assumptions, we can model the 
Bose-Hubbard system as a function of the lattice depth V as 

N N N 

H S iv) = j>icn + ~ E u *<y) - E Uw sVi+ Jjiv) &]_&), (i) 

j =1 j =1 j =2 

where ujj is the external site energy on j, Oj (dj) creates (annihilates) a particle on site j, Uj is the on-site interaction 
energy and Jj is the nearest neighbor tunneling matrix element between sites j — 1 and j. 

The deep lattice assumption allows us to approximate the dependence of the system parameters on V by treating 
the wells as coupled harmonic oscillators |29| . Under these assumptions, and in a frame rotating with the ground 
state energy of the perpendicular trapping potential, the Bose-Hubbard parameters have the following dependence on 
V : 


= / dx ^^ ( 2m + ^ lat ^ 

Uj(V) = ^t f d s? *)| 4 ( 2 ) 

m j 

J j(V ) = dx V$(aO + V lat {x)^j 


where ipj(x) is the ground state harmonic oscillator solution of the local lattice potential centered at site j, a s and 
m are the scattering length and mass of the atom, and 'ip±(y,z) is the harmonic oscillator solution of the trapping 
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potential. Analytic solutions for uij(V) and Uj(V ) are 

«>j(V) = 2 y/ErVv + V ex t(xj ) 

ETj(V’) = v^—£> /4 \/V> 1/4 (3) 

a l 

where all of the parameters are expressed in terms of recoil Energies E r = h 2 tt 2 / (2maf) . Throughout this work we 
assume a lattice spacing of a/ = 640 nm, the mass for 87 Rb whose scattering length is taken to be a s = 5.2 nm, and 
the perpendicular lattice depth is V± = 50E r . In our model, we obtain Jj(V) numerically. 

Previous works proposed by Refs. p~8l [20j have demonstrated that sinusoidal modulation of V(t) can lead to higher 
Bloch band excitations. However, assuming that V(t) ~ 15 E r guarantees a separation of parameters uj V J. 
Since the modulations presented below will consist of frequencies a < U, higher Bloch band dynamics can be safely 
neglected. This has been verified in our numerical simulations of higher band dynamics these systems. 


B. Time-Dependent Born-Markov Master Equation 


We study the transport properties of these one-dimensional systems using the tools developed in [50], where the 
ends of the lattices are coupled to reservoirs acting as sources and sinks of ultracold atoms as in Fig. [l] We use a 
quantum master equation approach pm mg to model the dynamics of the Bose-Hubbard system coupled to reservoirs 
of ultracold neutral atoms. Working in the zero-temperature limit with strongly repulsive atoms, the reservoir’s level 
occupancy can be characterized by a chemical potential /i, such that all states below /i are occupied and all states 
above are vacant. A chemical potential difference between the left and right reservoirs [Al > /iR, may induce an atomic 
current across the lattice from left to right. When U J the particle manifolds of the system become well-separated. 
If the values of Pl,r do not resonantly overlap with any system eigenenergies, the care taken in Ref. m is not 
necessary, and we can employ the standard Born and Markov approximations [321434] . Under these conditions, the 
value of fiL determines the maximum allowed particles on any given site. In our model, we assert the value of /al by 
truncating the system’s basis at the appropriate level. 

Since the modulation frequencies a will be of the same order of magnitude as the other frequencies associated with 
the system, the rapid decay of the memory kernel ensures that the time-dependence of the system has a negligible 
affect on the Liouvilian. This generates a master equation of the standard Lindblad form m- For an N-site lattice 
where atoms are being pumped onto site 1 by a reservoir with chemical potential /jll, and removed from site TV by a 
reservoir with chemical potential /ar = 0, our master equation model is 


dp 

dt 




^ (aia\p + paia\ - 2a\paij - | [a} N a N p + pa ] N a N - 2a N pa t N 'j , 


(4) 


where Hs(t ) = Hs(V(t )), and r is the decay rate of the reservoirs, which we assume to be identical for convenience. 


C. Transport Current 

We define the atomic transport, or current, in these systems by monitoring the average number of particles leaving 
the system. This can be determined from the rate contributions of the master equation. The reservoir coupled to 
the N th lattice site is responsible for particles leaving the system. This reservoir action is quantified in the last line 
of Eq.Q. In Fig. [2j we illustrate how the action of each of these Liouvillian terms affects the system populations. 
Taking the inner product of Eq. with a system eigenstate | j) who belongs to the n particle manifold of the system 
(for example) shows that the K(j\aNpd^ N \j) term corresponds to positive contributions, populating the state | j) from 
the n + 1 particle manifold, while the — K/2(j\a f N aNp\j) and —K,/2(j\cL [ N a]srp\j) terms lead to depletion of state | j) to 
the n — 1 particle manifold. Therefore, the net current can be defined as either the total populating of the system 
from all states above, or the total depletion to all states below. We choose the former, yielding: 

(7) = ^ K{j\a N palf\j) = /tTrfajvpat] = K(atajv) , (5) 

j 

where we have exploited the cyclic property of the trace, and recognize 

I = K {a^ N a N ) ( 6 ) 

as the current out of the system. This definition makes sense, since k is the decay rate of the Lindblad master 
equation. Thus the current out should equal the average number of particles on the N th site, multiplied by the rate 
at which the particles leave the system. 
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FIG. 2. An illustration of the actions of the reservoir connected to the N th lattice site on an n-particle system state. The 
CLNpCL jy term contributes to population of the n-particle state from the (n + 1) particle manifold, while the a^aNp and pa^ciN 
terms deplete population to the (n — 1) particle manifold. 


III. SINGLE ATOM EXCITATIONS 

We first consider cases where U J, and /il is set to maintain an occupancy of one atom on the first site. 


A. Modulating the Optical Lattice 

Consider a stationary (V(t) = V m [ n ) N-site system defined by Eq. Q with arbitrary external energies Uj. When 
energy differences between at least one pair of adjacent sites are large compared to their tunneling amplitudes, 
transport is greatly suppressed. Rotating this system with 

N 

H 0 = ^2tv j a]d j (7) 

3 = 1 


yields 


m = 5 


N 


Ujd]d)djdj 

3 = 1 


N 

E 

3 =2 


{^3 e 


iSnt 


d^dj 


-i + e 


jt 


)■ 


(8) 


where Sj = ujj — Uj-i, is the energy difference between sites j and j — 1. 

If each tunneling rate in Eq. © could be modulated with 

Jj Jje lSjt , (9) 

then the offset modulated system of Eq. Q is equivalent to a flat stationary lattice system, which would yeild the 
optimal current response. The modulation in Eq. © is ideal, since it would fully recover the maximal transport re¬ 
sponse from an arbitrary offset stationary lattice. The lesson learned from this example is that sinusoidal modulations 
of the tunneling rates may greatly enhance transport in lattices with energetic mismatches. 


B. Optimization Scheme 

Given an arbitrary 7V-site lattice, we propose that if there are M < N unique, nonzero values of \Sj\, that a great 
deal of the intrasystem current can be recovered if the overall lattice depth is varied such that the tunneling barriers 
experience a sinusoidal modulation consisting of a superposition of the M unique frequencies. 
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For example, 


V(t) = E min - - In 


1 


M 

— cos 2 

M ^ 

k=0 




( 10 ) 


would yield 


N 

Hs(t) *5>(t) 
./'=! 


1 JV 1 M 

2 E a } a } a j a i ~mJ2 cos2 

j=l &=1 


I* 


N 

E + Hx - 

i =2 


( 11 ) 


Here, f3 is the parameter from an exponential fit of the tunneling J(E) ~ Jm a xe _ ^ v_ymin ). Assuming H m i n = 15E r , 
we find J max — 6 x lO -3 ^ and /3 cx 0.24/E r . In order to compare results against the response of a stationary 
ideal lattice with a tunneling rate of J max , we’ve included the 1/M term to guarantee that the tunneling rates in the 
modulated Hamiltonian do not exceed J max . 

Although a clever choice of the relative phases of the oscillatory drive can enhance the current response even further, 
we adopt cosines in this work to yield a simpler analytic form. 


C. Effective Stationary Models 


In the deep lattice regime, stationary effective flat lattice models can be generated that accurately reproduce the 
dynamics of Eq 0 These effective models are significantly less demanding to numerically evolve. They also shed 
light on the transport recovery mechanism proposed above. 

Assuming the global modulation from Eq. 


(10), expressing the onsite interaction energy as a modulation about its 
mean U(t) = (U) + A U(t) and transforming into an interaction picture that rotates with the time-dependent version 
of Eq. ([7]) yields 


1 

Hs(t) =- £ «[/,) + A Uj(t)) 
3 = 1 



N 

x J jW Sjt & ] & 3- 1 + H - c - 
3= 2 


( 12 ) 


where 5j(t) ~ Sj in the assumed parameter regime. 

When two adjacent sites have the same energy, the stationary portion of the tunneling matrix element between 
these sites is J/2 . When Sj ^ 0, the stationary portion becomes J/(4M). Since \Sj\ A[/, J a secular approximation 
similar to the rotating wave approximation can be made: for timescales large compared to %/ min{|^ |}, the oscillating 
terms in Eq. |~L2| ) become insignificant compared to their stationary counterparts. Neglecting these terms yields 


where 


1 N 

= — y^(Uj)OjOjajdj 

N 

- OjUj-i + H.c. 

3= 1 

3= 2 


J 

II 

0 

1 

Jj = - 
J 2 

5 3 ± ° => 

T J 

j ~ 4M * 


(13) 


(14) 


That is, the assumed modulation of the gapped system yields an effective stationary Hamiltonian whose flat energy 
landscape assures the optimal current response. Although the rescaling of the tunneling rates implies a current 
response less than the corresponding flat lattice with uniform J, the current gained in these systems is still several 
orders of magnitude greater than the unmodulated lattice with one or more energy offsets. 
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D. Numerical Results 


In this section, we numerically verify transport recovery of Eq. 0 due to the modulation scheme presented in 
Sec. |IIIB[ and the accuracy of the effective stationary model in Sec. ImUl for a variety of five-site lattices. We 
first consider the system depicted in Fig. [3ja). Working with a minimum lattice depth of V m { n = 15 E r , we assume 
Sj = [—0.1, 0.3, —0.4, 0.2]E r . This lattice has four unique energy offsets. When coupled to reservoirs whose decay 
rates satisfy Jj/n — 15, this lattice exhibits a steady-state current response that is reduced by a factor of 1.3 x 10 9 
compared to the maximum response determined by its flat counterpart. Since we are assuming that /hl is set to 
maintain an occupancy of one particle on the left site in this section, the values Uj are irrelevant since the onsite 
interaction terms are not accessed. 

The effective stationary model for the dynamics of this lattice under the proposed modulation is illustrated in 
Fig.§b). A numerical simulation of the full unapproximated modulated dynamics of Eq. ( f~L2| ) for this energy-offset 



/ 2 J 3 J A J 5 

v/WW 


FIG. 3. The mapping of the five-site optical lattice depicted in Fig. [I] (a) onto its effective stationary Hamilton when modulated 
(b). In this example Jj = Jj/16- 


lattice is provided in Fig. [4] The current response is normalized by the steady-state maximal response of the flat, 
stationary lattice with uniform J. In steady-state, modulation leads to a current gain of 9.8 x 10 8 compared to the 
stationary offset lattice, which is 78% of the current response of the ideal system. The effective stationary model of 
Eq (13) is also plotted in this figure. For times large compared to timescales set by the energy differences of this 
system (t ~ Ti/E r in the plot), the unique dynamics of the modulated system are accurately recovered by the effective 
model, yielding a percent error of 0.2% from the full steady-state dynamics. 



t Er/fl 

FIG. 4. Time evolution of the current response of the modulated, offset lattice depicted in Fig. [3^a) (red solid) and the 
stationary effective Hamiltonian in Fig.plb) (blue dashed). The modulation increases the lattice’s current response by a factor 
greater than 9.8 x 10 8 , recovering over 60% of the lattice’s maximal current response. The stationary effective model accurately 
captures the evolution for times large compared to t — K/E r . 

The steady-state current gains and accuracy of the corresponding effective models for a variety of five-site lattices 
are provided in Table [IJ The characteristics of the Hamiltonian described above are presented in the fourth row of 
this table. In all four cases, the atomic current is largely recovered and the effective models accurately predict the 
dynamics of each modulated system. 

As E m i n gets reduced, the proposed modulation scheme is still capable of recovering most of the current. For 
example, assuming that V m [ n = 5 Er, J/k ~ 29 and that a Feshbach resonance can increase the scattering length 
















Sj(E r ) 

Current Gain 

Percent Recovered 

Percent Error of i7 e ff 

[0,0.1,0,0] 

5.6 x 10® 

64% 

1 . 0 % 

[0.2,0, -0.2,0] 

3.3 x 10 7 

46% 

1 . 5 % 

[0,0.1,0.3,-0.3] 

4.7 x 10® 

13% 

0 . 1 % 

[-0.1,0.3,-0.4,0.2] 

9.8 x 10® 

78% 

0 . 2 % 


TABLE I. The current characteristics of a variety of five-site lattices. The columns display the lattice characterized by its 
energy gaps, the steady-state current gain from modulation, the percent recovered compared to an ideal system and percent 
error of the corresponding effective Hamiltonian. 


so that the Uj terms remain irrelevant, simulations show that the Sj = [—0.1, 0.3, — 0.4, 0.2]E r lattice experiences a 
current gain of 4.7 x 10 5 where over 71% of the maximum current is recovered. The percent error of our effective 
model in this case is 30%. The effective model is breaking down since we are starting to push the limits of our 
approximation, but the current gained in the system remains significant. 


E. Gain as a function of maximum lattice depth 

In an experiment, it might not be advantageous or possible to vary the tunneling from some maximal J to zero. 
Next we examine the current response of the Hamiltonian depicted in Fig. |3ja) as a function of the maximum lattice 
depth V max assuming V m i n = 15E r . Figure [5] presents this simulation where the current response is normalized to its 
steady-state value when Vm ax = 50 E r . The data from this plot shows that varying the lattice depth by 2 E r , 5 E r , or 



FIG. 5. Current response of the modulated lattice in Fig. |3|a) as a function of Vmax assuming V m in = 15 E r . 

SE r recovers respectively 38%, 79%, or 91% of the maximum current obtained when V max = 50 E r . Thus, even small 
depth changes lead to significant current gains. 


IV. MULTIPLE ATOM EXCITATIONS 

The optimization condition changes when the left reservoir is set to maintain an occupancy greater than one particle 
on the left site. We close by considering a two-site lattice with a finite energy offset S 2 , where /il is set to maintain 
an occupancy of two atoms on the first site. For an arbitrary modulation frequency <a, the system Hamiltonian we 
consider is 


H(t) — + ^i(t) + £2 

(U) + A U(t) 


a\a 2 


djdjdidi 


■ a\a\a 2 a 2 


— J cos 2 ( — 


(!*) ( 


a\ai 


a\a 2 


)■ 


(15) 
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Assuming F min = 15 E r and F max = 50 E r implies S 2 (t) ~ S 2 , (U) \AU(t)\ and ( U) / J ~ 390. We futher assume 

S 2 = 0.17% and J j 15. 

The numerical steady-state current response of Eq. ( p~5| ) in Eq. <@ as a function of the modulation frequency a 
is presented in Fig. [6j In the previous section, the current was greatly enhanced when the system was modulated 



FIG. 6. Numerical current response for an energy mismatched two-site lattice driven to maintain an occupancy of two atoms on 
left site. Modulating at S 2 (corresponding to ~ 0.2 on the horizontal axis), leads to a current response that is small compared 
to the one obtained when modulating the lattice at S 2 — ( U ). 


at the frequency of the gap S 2 . Here, when a = S 2 , the current increases by a factor of 7.3 x 10 3 . However, when 
a = S 2 — (U ), a current gain of 1.7 x 10 5 is observed. The effective stationary Hamiltonian corresponding to Eq. (15) 
with a = S 2 — (U) is 


#eff = cJidjdi + (ui + {U))a\a 2 + + alala 2 a 2 ^j - ^ (a 2 ai 


+ Cl\cL2 


■) 


(16) 


The percent error between current predicted by this effective model and the full modulated system is 7.2%. 

The change in the optimization condition can be understood by examining how the energetic relationships of the 
unmodulated system’s Fock states effectively change when the system is modulated at frequencies S 2 and S 2 — (U ). 
The analysis that follows is closely related to that which was originally done in Ref. [3]. Figure [7] depicts the bare Fock 
energies of the unmodulated, energy-offset lattice. In the absence of any modulation, there are no system degeneracies, 
and the action of the reservoirs evolve an arbitrary system almost completely into the |2,0) state. 


12 , 2 ) 


-[ 


| 1 , 2 > 


12,1 > 


S 2 + {U) 
S 2 - (U) 




| 0 , 2 ) 


|2,0> 


|U> 




| 0 , 1 > 


IhO) 


| 0 , 0 ) 

FIG. 7. An illustration of the relative relationship of the uncoupled Fock states, for zero to two particles per site, of a two-site 
lattice with finite 5 2 and (U). 
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Figure |8ja) displays the energies for the effective Hamiltonian obtained by modulating the system at 82 - This 
establishes a resonance between 11, 0) and |0,1), as well as |2,1) and |1,2). When restricted to the single particle 
excitation case as in Sec. m this condition optimizes transport since the action of the reservoirs drives population 
cycles through the lowest four states in the diagram. However, when /il is set to maintain an occupancy of two 
particles, the energy mismatch between |2,0) and 11,1) assures that, regardless of the system’s initial condition, the 
action of the reservoirs evolve the vast majority of the system population into the 12, 0) state. Current may leave 
the system, but only via a second-order, off-resonant transition from |2,0) —» |0,2). Such a transition leads to the 
relatively small current gain seen in Fig. <©> when a ~ 0.2 x (*2 - (U)). 



FIG. 8. Relative energies for H e ff of a two-site lattice with finite £1 and U when it is modulated at (a) 82 and (b) 82 — ( U }. 
The blue arrows represent actions of the left reservoir, the red arrows represent actions of the right reservoir, while the green 
arrows represent intrasystem transitions. As seen in (a), modulations with 82 create an effective resonance between 11,0), and 
10,1). However, the energetic mismatch between |2, 0) and 11,1) prohibits effective transport when the reservoir is driving the 
|2,0) state, (b) Alternatively, if the system is modulated with 82 — (U), the resonance between |2,0) and 11,1) recovers the 
maximal transport response for this system. 


Figure |8jb) shows the energetic relationship of Eq. (16), the effective Hamiltonian associated with modulating the 
system at 82 — (U). This modulation makes the system in question equivalent to the forward-biased atomtronic diode 
introduced in Ref. 0, but with J reduced by a factor of 4. The resonance between the |2,0) and 11,1) states assures 
a large transport response via cycles between the 11,0), 12, 0), |1,1) and |2,1) states. 

Thus appreciable transport enhancement in higher particle manifolds is possible, but care needs to be taken to 
ensure that the correct resonance conditions are met. 


V. CONCLUSION 


In this article, we have demonstrated that varying the depth of an optical lattice so that the tunneling barriers 
experience a sinusoidal modulation can greatly enhance the transport characteristics that are otherwise suppressed 
by energetic mismatches in the system. Through a secular approximation, we have shown that this modulation maps 
a stationary lattice that prohibits transport onto an effective stationary lattice that allows for optimal transport. 
We have demonstrated transport recovery for both single- and double-site occupied lattice sites. For the case where 
at most one atom was allowed on any given site, we verified the behavior of our proposed modulation schemes as 
well as our effective stationary models for a variety of five-site lattices. We also examined the current response as a 
function of the modulation intensity and concluded that even very small variations of the lattice depth can generate 
significant gains in transport. Finally, we demonstrated how current can be recovered when higher particle manifolds 
are accessed. For this case, the conditions for optimizing transport change: the modulation frequencies need to be 
chosen to create the correct effective resonances across particular particle manifolds in the system. 
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